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Abstract: The discovery of a Standard Model Higgs boson is possible when experimental 
cuts are applied which increase the ratio of signal and background cross-sections. In this 
paper we study the pp — > H — > WW signal cross-section at the LHC which requires a 
selection of Higgs bosons with small transverse momentum. We compare predictions for 
the efficiency of the experimental cuts from a NNLO QCD calculation, a calculation of the 
resummation of logarithms in the transverse momentum of the Higgs boson at NNLL, and 
the event generator MC@NLO. We also investigate the impact of various jet- algorithms, 
the underlying event and hadronization on the signal cross-section. 
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1. Introduction 



Physical processes at collider experiments can be simulated using flexible event generators 
such as PYTHIA [1] and HERWIG [2]. In this approach, the momenta of partons in hard 
scattering processes are distributed among hadrons using an approximate probabilistic 
algorithm for parton branching and hadronization. Traditionally, event generators compute 
the hard scattering partonic cross-section at leading order in fixed order perturbation theory 
which yields only a rough estimate. 

Cross-sections for the hard interaction of three or four particles can be computed 
routinely through next-to-leading order (NLO) in the ag expansion. The results from 
NLO QCD calculations and parton shower event generators are often combined with an 
empiric method. First the efficiency of the experimental cuts and normalized differential 
distributions arc computed with parton shower event generators. Then, they are multiplied 
with the result for the total cross-section from the NLO calculation ^. However, it is now 
understood how to combine parton shower generators and NLO results with theoretically 
sound methods so that (i) leading logarithms are resummed with the parton shower, and 
(ii) all differential cross-sections are exactly accurate through NLO upon an expansion of 
the Sudakov factors in [5,6]. 

Fixed order perturbativc computations have now advanced beyond the next to leading 
order, and there are two hadron collider processes which are known through next-to-next-to- 
leading order (NNLO) in ag. The NNLO total cross-section of the Drell-Yan process [7,8] is 
the most precise theoretical prediction for a hadron collider observable with a scale variation 
uncertainty of about 1%. The total cross-section for Higgs boson production is also known 
at NNLO [8-12]. For the LHC, the NLO [13,14] and NNLO corrections are both important 
and increase the LO result by about 70% and 30% respectively. The perturbative series 
converges slowly and the remaining theoretical uncertainty is about ±10% [8-12,15]. 

In Higgs boson [16-18] and electroweak gauge boson [19-22] production there exist 
novel differential cross-section calculations at NNLO [16-22] ^. The cross-sections with 
arbitrary experimental cuts applied at the parton level can be computed exactly at this 
order in pertubration theory for the two processes. It is very instructive to compare the 
efficiencies of experimental cuts from the newly available fully differential NNLO calcu- 
lations, merged NLO and parton-shower calculations with MC@NLO, and simple leading 
order event generators. This is valuable in order to estimate the inherent theoretical un- 
certainties of the above approaches. 

Such a comparison can be made for the Drell-Yan process with the results from 
Refs. [21,22,26], as well as discussed in Ref. [21]. MC@NLO and the NNLO Monte-Carlo 
FEWZ [21,22] predict very similar experimental efficiencies for the entire kinematic range 
where the NNLO prediction retains its phenomenal scale variation of about 1%. Signif- 

^Note that differential factors" to reproduce bin integrated differential distributions at higlier orders 
in perturbation theory are also used [3,4]. 

^In non-hadronic collisions the state of the art at NNLO is fully differential cross-sections for e^e~ — > 
3 jets [23-25]. 
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icant differences, however, arise when the experimental cuts suppress contributions from 
the two-loop matrix elements. 

A similar comparison [27] between MC@NLO and the NNLO partonic Monte-Carlo 
FEHiP has been made for the Higgs boson diphoton signal at the LHC. The signal cross- 
section is known with a scale uncertainty of about ±7% at NNLO [17]. Besides a relatively 
large perturbativc correction from NLO to NNLO of about 20%, the efficiency of the 
experimental cuts turns out to be very similar in MCONLO and NNLO [27]. 

A challenging channel for the Higgs boson discovery at the LHC is pp — > H — > WW. 
This channel is contaminated by background processes, pp it and pp WW, with 
much larger cross-sections. For a Higgs boson with a mass close to the W-pair thresh- 
old, the cross-sections for all other discovery signals are suppressed. In this case, the 
pp ^ H ^ WW lulu process becomes the only viable channel for the Higgs boson to be 
discovered at the LHC. It is thus indispensable to achieve a very good signal to background 
(S/B) ratio. An optimized selection of W-pair events [28] is then required. If appropriate 
cuts are applied, as in Refs. [29-32], a discovery of a Standard Model Higgs boson with a 
mass close to the threshold should be achieved with an integrated luminosity of a few fb~^ 
at the LHC. 

The pp ^ H — > WW Iviv decay mode was recently implemented in the NNLO 
Montc-Carlo FEHiP [17] and a calculation of the cross-section with the experimental cuts 
of Ref. [29,30] was performed in Ref. [33]. The very good agreement of the MC@NLO and 
NNLO calculations for the efficiency of the experimental cuts in the diphoton signal [5,17, 
27] does not guarantee that a similarly good agreement will be found in the H WW chan- 
nel. In this paper, we will compare the NNLO predictions from Ref. [33] with resummation 
calculations and event generators. 

In the H ^ 77 channel, the experimental cuts select events "democratically" , irrespec- 
tively of the transverse momentum of the jets associated with the Higgs boson production. 
The event selection in the H WW channel [30, 31] imposes an explicit jet-veto and 
other cuts which reject events with large transverse momentum pi^ for the Higgs boson. 
This may turn out to be problematic for an agreement in experimental efficiencies between 
MC@NLO and NNLO for two reasons. First, the NNLO/NLO i^-factor is sensitive to the 
cutoffs imposed on the [16,17,34]; this effect is treated only in the parton shower approx- 
imation in MC@NLO. Second, by selecting events with low ]9^, multiple gluon radiation 
effects which are not included in the fixed order NNLO calculation may be important. 

The resummation through next-to-next-to-leading logarithms (NNLL) of pi^ is now 
achieved [35-37]. The resummed spectrum, after matching to fixed order, integrates to 
the NNLO total cross-section. This calculation takes into account consistently both mul- 
tiple gluon radiation effects at low transverse momentum and fixed order high transverse 
momentum contributions. 

A theoretical prediction for the signal cross-section pp — >^ H — >^ WW ivtv cannot be 
made directly from the spectrum, since the experimental cuts restrict many phase-space 
variables. The cross-section must therefore be computed using fully differential Monte- 
Carlo programs. However, we will use the resummation calculation of the pJjJ spectrum [36- 
38] to validate fixed order Monte-Carlo's and parton shower event generators in the low 
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kinematic region, which is favored by the experimental selection cuts. 

We first compare in Section 2 the theoretical calculations from the NNLL resumma- 
tion [38], MC@NLO [5] and NNLO [33] for the spectrum. Then we compare in Section 3 
the MCONLO [5], HERWIG [2] and NNLO [33] Monte-Carlo predictions for kinematic dis- 
tributions of variables which are restricted by the experimental selection.^ We perform a 
similar comparison for the signal cross-sections when all cuts are applied. We find a good 
agreement for the efficiencies of experimental cuts and normalized kinematic distributions. 
This gives us confidence that the selection of events for the signal cross-section does not 
invalidate the approximations in the used Monte-Carlo programs. 

In Section 4 we study the sensitivity of the signal cross-section with all experimental 
cuts applied to the choice of jet algorithm. We find a mild change on the cross-section 
(~ 6%) by using the SISCone or the fcT-algorithm. We also study the effect of hadroniza- 
tion (using the model in HERWIG [2]) and of the underlying event (using the JIMMY 
model [39]). We find that proposed experimental cuts render the cross-section sensitive 
to these effects at the level of up to 5 — 10%. The changes to the cross-section are with 
opposite signs and the combined effect is rather mild. 

2. Integrating over the transverse momentum of the Higgs boson 
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Figure 1: Cumulative cross-section for the Higgs transverse momentum distribution at fixed NLO 
and NNLO expansion. 



In this Section we compare fixed order [17] and parton shower calculations [2,5] against 
the theoretical predictions for the transverse momentum spectrum of the Higgs boson 

^In the case of MC@NLO and HERWIG, we use a modified version of HERWIG with the correct 
decay angular correlations. 



-3- 



from resummation [38]. Fixed order perturbation theory is invalid for small values of 
p^. Nevertheless, observables can be reliably computed if their definition allows for an 
integration over a sufficiently large range in values. Here, we study the cumulative 
cross-section 

H.max 

Pt 

a{p^'—)= I ^dp^. (2.1) 


This observable mimics the effect of selection cuts with a cutoff on the maximum Higgs 
transverse momentum. Since in the transverse plane the Higgs boson balances the asso- 
ciated jet radiation, the cross-section with a jet-veto is similar to the cross-section with a 
veto on high values of p^. 



pp ^ H + X 
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Figure 2: Cumulative cross-section for the Higgs transverse momentum distribution with 
MC@NLO, NLL resummation [38] and at NLO. 



In Fig. 1 we plot the cross-section in the fixed NLO and NNLO expansion. As in all 
plots of this paper, we show the scale variation in the interval m}i/2 < = fiji < 2mH. 
The fixed order cross-sections at NLO and NNLO become negative and tend to infinity 
when the p^ cutoff is small (below 10 GeV ). For such small values of p^'™^^ perturbation 
theory breaks down. We observe that for larger cutoffs (above 40 GeV and up to 100 GeV) 
the NLO and NNLO results are in very good agreement. The NNLO cross-section increases 
however faster than NLO with higher cutoffs leading to the known by ^ 20% larger NNLO 
result with respect to NLO for the total cross-section [8-10]. The analyses in Ref. [29,30] 
show that a better S/B ratio is achieved if stricter than 40 GeV cutoffs are used for the 
jet transverse momenta; for cutoff values in between 20 GeV and 40 GeV we observe large 
perturbative corrections. 
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In Fig. 2 we compare the fixed NLO result with the resummed NLL spectrum from 
Ref. [38] and with MC@NLO [5] without hadronization and underlying event. We find that 
for < 40 GeV the parton shower or the NLL resummation change significantly the inte- 
grated NLO distribution. All results converge for higher and higher cutoffs pJ^'™^^ and 
agree with each other for the fully inclusive cross-section. Notably, the MC@NLO and the 
NLL resummation are in a rather good agreement with each other within the uncertainties 
from scale variation. 

Differences between cross-section predictions with fixed order perturbation theory and 
resummation are expected to become smaller when the fixed order calculations are extended 
to higher orders. In [33] it was found that the average of the Higgs boson, when 
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Figure 3: Cumulative cross-section for the Higgs transverse momentum distribution at NNLO in 
fixed order and with NNLL resummation [38] . The two approaches agree very well in the kinematic 
range which is relevant for the envisaged experimental cuts. 



discovery selection cuts are applied, can be as low as < >~ 15 GeV. This leads to 
the question whether the NNLO result, unlike the NLO result, is a reliable prediction for 
such small values of the average p^. In Fig. 3 we compare the integrated p^ distribution 
at NNLO against the resummed NNLL spectrum. We find a very good agreement between 
the two approaches for surprisingly low values of p^'™''^. We conclude that higher than 
NNLO perturbative contributions, which are accounted for with the resummation, remain 
small for the integrated p^ spectrum, when the maximum Higgs transverse momentum is 
restricted even down to 20 GeV. 

We have now validated the NNLO perturbative calculation [33] in a challenging case 
for fixed order perturbation theory which similarly emerges, due to the jet-veto and other 
cuts favoring small p^ values, in the search for a Higgs boson in the WW decay channel. 
However, the simulation of processes at NNLO is only performed at the parton level. We 
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would like to investigate whether parton shower Monte-Carlo programs, which can also 
model non-perturbative effects and are computationally more ffexible than NNLO Monte- 
Carlo's, provide realistic estimates of the signal cross-section. 

We first discuss the problem of the normalization of the event generators. Parton 
shower Monte-Carlo programs predict the same total cross-section as the cross-section for 
their encoded partonic hard scattering at fixed order in perturbation theory. Therefore, 
HERWIG predicts the Higgs boson total cross-section with LO accuracy (underestimating 
it by a factor of ~ 2) and MC@NLO provides NLO precision (underestimating the total 
cross-section by a factor of ~ 1.25). A matching of parton showers to NNLO fixed order 
calculations is not yet developed. Following a practical approach, we will validate whether 
the efficiency of experimental cuts and normalized differential distributions are in agreement 
with the NNLO calculations of Ref. [33]. We will then rescale the predictions of the 
MC@NLO and HERWIG event generators with a global ii'-factor in order to reproduce 
the fixed order result for the total cross-section. We will denote that the results of the 
Monte-Carlo X have been multiplied with a iC-factor using the notation R(X). 



pp^H+X pp^H+X 




Figure 4: Cumulative cross-section for the Higgs transverse momentum distribution. The scaled 
MC@NLO and HERWIG spectra agree very well with the resummed NNLL spectrum [38]. 

Now we will test how well event generators agree with resummation results for the 
spectrum. In Fig. 4 we compare the integrated spectrum of MC@NLO and 
HERWIG against the resummed NNLL prediction. We observe that both generators 
are in very good agreement with the NNLL spectrum. This is especially surprising for 
HERWIG which aims to describe the salient physics features of the process. Note, how- 
ever, that MC@NLO gives slightly larger and HERWIG slightly smaller values than the 
NNLL resummation [38]. 

Before we conclude our analysis of the integrated distribution we wish to comment 
further on the scale variation of the fixed order results. In Fig. 1 we find a pi^'™'^^ with 
no scale variation. A similar behavior is also observed for the accepted cross-section with 
all experimental cuts [29,30] in Ref. [33]. It is therefore necessary to investigate in better 
detail the variation of the NNLO result. We show the value of the cross-section at NNLO 
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for prp'™*^ ~ 37GeV, varying independently the renormalization and factorization scales 
(the errors correspond to the numerical integration): 



a \ph\ 




f^F = ^ 

HP = mn 
HF = "^rnu 


13.31 ±0.13 13.76 ±0.08 13.45 ± 0.05 12.82 ± 0.04 
13.15 ±0.13 13.85 ±0.08 13.69 ± 0.06 13.14 ±0.04 
13.13 ±0.13 14.00 ±0.08 13.96 ± 0.06 13.47 ± 0.04 
13.05 ±0.13 14.15 ±0.08 14.21 ± 0.06 13.76 ± 0.04 



The cross-section is more sensitive to independent variations of the renormalization and 
factorization scales; however, this variation is significantly smaller than in the total cross- 
section. A detailed study of the cross-section when all experimental discovery cuts are 
applied, which shows a similar scale- variation pattern, has been made in Ref. [33]. 

Summarizing, in this Section we found that the integrated spectrum is predicted 
reliably at NNLO for the kinematic range of ]^ which is relevant in the search pp — > H ^ 
WW iulv. On the contrary, the NLO fixed order calculation is unreliable. We have 
also established that for the same observable, MC@NLO and HERWIG are in a very good 
agreement with the resummed NNLL spectrum when they are normalized to a common 
NNLO total cross-section. 

3. Kinematic distributions and signal cross-section 

The main backgrounds for the pp ^ H ^ WW — > Ivlu process are pp — > tt and 
pp WW. These backgrounds are sufficiently suppressed to allow for the discovery 
of a Higgs boson with a combination of experimental cuts [29,30], exploiting the spin- 
correlations in the decay of the Higgs boson and the high average transverse momentum of 
jets in top-pair events. In Ref. [33], the signal cross-section with these cuts was computed at 
NNLO. In this section we will compare the NNLO results of Ref. [33] with MCONLO. The 
public version of MC@NLO includes only partial spin-correlations for the H WW Ivlv 
decay. Here, the full spin-correlations for the decay of the Higgs boson have been imple- 
mented in MC@NLO. All the results of this paper correspond to a Standard Model Higgs 
boson mass of mn = 165 GeV. 

In Figs. 5-8 we present the cross-sections when a single cut is applied on 

• vau^ the invariant mass of the charged lepton pair, 

• 4>U,^ the angle between the two charged leptons in the plane transverse to the beam 
axis, 

• ^'Tmax' transverse momentum of the harder lepton, and 

• E^^^^, the missing transverse energy. 

The four distributions show an excellent agreement for the efficiencies among NNLO and 
MC@NLO. This is a remarkable result and could not have been easily foreseen; there is 
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pp^H+X^WW + X^e^yei/ + X 
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mr [GeV] 

Figure 5: Cross-section wlien tlie lepton invariant mass is constrained in tlic interval [12 GeV, w^"*] 
at NNLO and with MC@NLO. 



pp ^ H + X ^ WW + X ^ e""!/ e"i/ + X 




Figure 6: Cross-section for the transverse opening angle of the two leptons in the interval [0, (/i^^']. 



a significant change in the shape of the (pu, Pxmax' ^^'^ E^'^^^ distributions from NLO to 
NNLO, as seen in Ref. [33]. 

A crucial experimental cut for suppressing the top-pair contribution to the background 
is a jet-veto. We veto events which have a transverse momentum of the leading jet in the 
central rapidity region (|?/jet| < 2.5) that is larger than p^^°. For the jet definition we 
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pp ^ H + X ^ WW + X ^ e""!/ e"i/ + X 
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Figure 7: Cross-section wlien tlie maximum transverse momentum of tire leptons is in the range 
30 GeV,pi^™*j^ . Each lepton should have a transverse momentum of at least 25 GeV. 
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Figure 8: Cross-section when the missing transverse energy is i?™'^^ > E^^^^"™^. 



use here a algorithm with a jet-radius parameter R = 0.4; later we will also use a cone 
algorithm (SISCone [40]). The two algorithms are identical for the LO and NLO calculation 
where only up to one parton can be present in the final state, if the same jet-radius R is 
used. They differ, however, in the parton shower calculations (MC@NLO and HERWIG) 
and at NNLO as more partons are generated. In the envisaged experimental analysis a 
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jet-veto with a rather smah value of '^ 25 — 40 GeV is considered. We wih investigate 
whether the NNLO and MC@NLO predictions are consistent with each other for such small 
values of the jet- veto. 

The good agreement of the integrated distribution between NNLO, MC@NLO and 
NNLL resummation suggests that a good agreement between MC@NLO and the NNLO 
cross-sections with a jet-veto may also hold. The jet-veto cross-section should be qualita- 
tively similar to the cross-section with a cutoff on the since at NLO the Higgs transverse 
momentum corresponds exactly to the transverse momentum of the additional jet. How- 
ever, the two cuts are not exactly the same and they compare only qualitatively. The 
jet-veto applies only at central rapidities; in addition, beyond NLO the is not the same 
variable as the maximum transverse momentum of the jets. In Fig. 9 we present the cross- 

pp ^ H + X ^ WW + X ^ e""!^ e"i/ + X 
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Figure 9: The Higgs production cross-section with a fixed-order computation (NNLO) and 
MC@NLO rescaled with an inclusive if-factor (R(MC@NLO)) when a veto on jets with pT > *° 
at central rapidities \ri\ < 2.5 is applied. 

section with a jet-veto applied. Indeed, we find a very good agreement between the NNLO 
result and MC@NLO (rescaled with the appropriate NNLO/NLO i^-factor for the total 
cross-section) . 

In Table 1 we list the cross-section after all signal cuts as described in Ref. [33] are 
applied. We have used both the /ct and SISCone algorithm of Ref. [40] and their imple- 
mentation from Ref. [41]. The jet radius in the azimuth-rapidity plane was set to -R = 0.4 
and the merging parameter for the SISCone algorithm to / = 0.5 ^. The two algorithms 
yield formally identical results for the fixed order calculation through NLO and indistin- 
guishable results at NNLO within our Monte-Carlo integration precision ^. In the first 

*The merging parameter / defines, how much two separate proto-jets need to overlap in order to be 
merged into one jet. 

^We thank Gavin Salam for pointing out to us that the SISCone and A;t algorithms are formally different 
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section of the Table we present the results obtained using a fixed-order LO computation 
and the LO+parton-shower event generator HERWIG. We find a much larger fixed order 
result than when a parton shower is added. At fixed leading order all events have a Higgs 
boson with zero transverse momentum, and the jet-veto rejects none of the events. On the 
contrary, HERWIG generates a large fraction of events with p!^^ > 



O-acc [fb] 

jet algorithm 


„ _ mn 

SISCone A;t 


/7, = 2 mu 
SISCone 


LO 

HERWIG 


21.00 ±0.02 
11.16 ±0.04 11.59 ±0.04 


14.53 ±0.01 
7.60 ± 0.03 7.89 ± 0.03 


NLO 

MCONLO 

(HERWIG) 


22.40 ± 0.06 
17.42 ± 0.08 18.42 ± 0.08 

19.79 ±0.07 20.56 ±0.07 


19.52 ± 0.05 

13.60 ± 0.06 14.39 ± 0.06 

14.61 ±0.05 15.17 ±0.05 


NNLO 

^NNLO (HERWIG) 


18.18 ±0.43 18.45 ±0.54 
19.33 ± 0.09 20.43 ± 0.09 
22.02 ± 0.08 22.88 ± 0.08 


18.76 ±0.31 19.01 ±0.27 
17.24 ±0.07 18.24 ±0.07 
18.65 ±0.07 19.38 ±0.07 



Table 1: Cross-sections after the signal cuts of Ref. [33] are applied for different calculation meth- 
ods. The statistical integration errors are shown expicitly. The MC@NLO and HERWIG cross- 
sections are evaluated with 1,000,000 generated events. The fixed-order results were computed in 
Ref. [33] and require the Monte-Carlo integration of multiple sectors [17]. 

In the second section of Table 1 we present the results obtained using a fixed-order NLO 
computation, the event generator MC@NLO and HERWIG after we have rescaled it with 
an inclusive NLO/LO factor. While in MCONLO and fixed order (LO, NLO and NNLO) 
we can set the renormalization and factorization scales for the hard scattering at will, we 
use HERWIG at the default scale since this affects the triggering of the hadronization 
procedure. We then rescale the HERWIG result using a K-factor, taking the NLO fixed 
order result to be at the scale fi = mH/2 or /x = 2mH. The NLO result is quite different from 
the one obtained with MC@NLO. We can attribute this failure of the NLO computation 
to the poor modeling of the low region, as is shown in Fig 2. 

In the last part of Table 1 we present the results obtained using a fixed-order NNLO 
computation and the results from MC@NLO and HERWIG, rescaled to the NNLO total 
cross-section. The NNLO result and the rescaled MC@NLO give consistent results, albeit 
with different behavior when varying the renormalization and factorization scales. A de- 
tailed analysis of the NNLO scale dependence when all cuts are applied can be found in 
Ref. [33]. 

Wc note that for all the results of this section we use MC@NLO and HERWIG at the 
parton level, switching off the hadronization and without using a model for the underly- 
ing event. We observe that the /cx-algorithm gives larger cross-sections than the SISCone 
algorithm for MC@NLO and HERWIG; as mentioned before, the results for the two algo- 

already at NNLO. 
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rithms are very similar at fixed order through NNLO (within our integration precision). 
Additionally MC@NLO gives slightly smaller values for the cross-sections than HERWIG. 

We have established in this section that the efficiency of experimental cuts computed 
with MC@NLO and HERWIG is similar to the efficiency obtained at NNLO. There are 
rather dramatic changes in differential distributions when going from NLO to NNLO [33] . 
It is only at NNLO that the fixed order calculation is consistent with the parton shower 
efficiency of the experimental cuts. In the following Section we will study the dependence of 
the cross-section on effects that are not captured by the fixed order NNLO calculation. We 
will study the effect of hadronization and of the underlying event. We will also investigate 
further the differences in the cross-section due to the two different jet algorithms. 

4. Jet algorithms, hadronization and the underlying event 

In this Section we perform a study of the signal cross-section with all cuts applied using 
MC@NLO. We will analyze the impact of different jet clustering methods, hadronization 
and the underlying event. 

In Fig. 10 we plot the cross-section using MC@NLO as a function of the value for 
the /ct and the SISCone algorithm with a jet-radius R = 0.4 and R = 0.7. The clustering 
is applied to all final state particles before hadronization. We find that for small values of 



pp ^ H + X ^ WW + X ^ e^i/ e"y + X 




10 20 30 40 50 60 

p;^'° [GeV] 

Figure 10: Comparison between the cone and fcx algorithm for different values of the allowed 
maximum jet transverse energy (jet- veto). All other cuts are set to the values chosen in Rcf. [33] 
for the signal cuts. 

the jet- veto parameter the choice of the jet clustering method is more significant. For a 
jet-veto at p^^° = 25 GeV the choice of jet-algorithm changes the cross-section by ~ 6% 
with MC@NLO. A similarly large variation of ~ 7% is observed when we vary the jet- 
radius from R = 0.4 to R = 0.7. For a jet veto value larger than about ~ 40 GeV 



- 12 - 



the sensitivity of the cross-section to the choice of the jet-algorithm or the jet-radius falls 
below ~ 2 - 3%. 




Figure 11: Difference of the cross-section after signal cuts including the underlying event and 
hadronization models, with respect to the partonic cross-section. The cross-section is shown as a 
function of the jet-veto value for the SISCone clustering algorithm. 



We now study the effect of hadronization as it is modeled in HERWIG and of the 
underlying event as implemented in JIMMY [39]. In Fig. 11 we present the relative differ- 
ence of the cross-section with respect to the partonic cross-section when the hadronization 
or/and the underlying event models are switched on. We have used here the SISCone 
algorithm with a merging parameter / = 0.5 and two values for the jet-radius R=0.4 (left) 
and R=0.7 (right). We apply the signal cuts set to the values which are used in Ref. [33]. 
We vary, however, the allowed maximum value of p^*. Of interest are values of the jet-veto 
between 25 and 40 GeV, which are envisaged in the Higgs boson search. 

Qualitatively, we anticipate that the hadronization and the underlying event change 
the partonic cross-section with opposite signs (we refer the reader to the recent analysis 
in Ref. [42] for a detailed study). Hadronization reduces the average px of (gluonic) jets 
by roughly 6pT ~ (lGeV)/i?. The underlying event increases the jet pT by roughly 
6pT ~ i?^ X (5 GeV) at the LHC. The slope of the partonic cross-section with the jet-veto 
cutoffs (Fig. 9) is large for small values of the jet- veto. The shifts 5pT from hadronization 
and the underlying event can therefore induce significant changes to the cross-section. A 
jet-veto after hadronization corresponds to a looser effective jet-veto at the parton level. 
We therefore anticipate the cross-section to increase by switching on the hadronization 
model. Similarly, we anticipate a decrease of the cross-section due to the underlying event. 

The trends can be verified in Fig. 11. A smaller jet-radius increases the impact of 
hadronization and decreases the impact of the underlying event. The two effects are not 
linearly additive. However, we find that a cancelation between the two effects, which varies 
according to the jet-radius, takes place. For a jet veto p^^° = 25 GeV and a radius R = 0.4, 
the hadronization shift is about ~ 7% and the underlying event shift is ~ 4%. For a larger 
radius R = 0.7, the two shifts are 5% and 10% correspondingly. 



-13- 



pp ^ H + X ^ WW + X ^ e*i^ e"i/ + X pp ^ H + X ^ WW + X ^ e*i^ e"i/ + X 




■ 10 ao 30 40 50 60 ' 10 20 30 40 50 60 



p'i"-' [GeV] p'i"-' [GeV] 

Figure 12: Difference of the cross-section after signal cuts including the underlying event and 
hadronization models, with respect to the partonic cross-section. The cross-section is shown as a 
function of the jet-veto value for the clustering algorithm. 

In Fig. 12 we show the effects of hadronization and the underlying event for the fcx al- 
gorithm. We find the same features qualitatively for the two effects as in the SISCone 
algorithm. We note however, that the /cx-algorithm shows an overall reduced sensitivity. 

5. Conclusions 

In this paper we have studied QCD effects for the process pp — > H ^ WW — > iuiv. This 
process is of particular interest at the LHC since it is possible that for a range of mass 
values of the Higgs boson, this channel is the only viable one for a discovery. 

The cross-section with all envisaged experimental cuts applied was computed in an 
earlier publication [33] at NNLO in QCD. In this paper we compared these NNLO results 
with the leading order event generator HERWIG [2] and the event generator MC@NLO [5] 
which performs a matching of HERWIG with NLO fixed order perturbation theory. We 
found very good agreement in efficiencies of all experimental cuts that are relevant in the 
search for the Higgs boson. This is rather spectacular given that there are significant 
corrections in the total cross-section and the shape of kinematic distributions from NLO 
to NNLO. 

The experimental cuts select events with small transverse momentum of the Higgs 
boson. This is important in order to reduce the selection of events from top-quark pro- 
duction, which is a major background. We have compared a NNLO computation and the 
result of NNLL resummation [38] for the cumulative distribution. We found that NNLL 
resummation does not induce significant corrections with respect to the NNLO calculation 
for the kinematic range which is favoured by the selection cuts. We have also found that, 
within the uncertainty from scale variations, MCdNLO and HERWIG are in very good 
agreement with the NNLL result. On the contrary, fixed order NLO perturbation theory 
provides a rather poor approximation for the required distributions and efficiencies. 
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Finally, we investigated the magnitude of effects due to hadronization, the underlying 
event, and the choice of jet algorithms. For typical choices of parameters and cutoffs in the 
experimental cuts we find a mild dependence of the cross-section on these effects. In this 
paper we did not examine a variety of remaining uncertainties, such as uncertainties in the 
parton densities, which may be relevant at a precision level of 10%. Studies of electroweak 
corrections can be found in Refs. [43-45] and of the background-signal interference in 
Ref. [46]. 

The results and comparisons made in this paper provide a firm validation of the fixed 
order NNLO results and the event generator tools which are available for simulating the 
pp ^ H — > WW — > luiv process at the LHC. 
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